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Abstract: We use a ten dimensional dual string background to study aspects of the physics 
of large N four dimensional SU{N) gauge theory, where its fundamental quarks are charged 
under a background electric field. The theory is M = 2 supersymmetric for vanishing tem- 
perature and electric field. We work in a limit where the quarks do not back-react. At zero 
temperature, we observe that the electric field induces a phase transition associated with the 
dissociation of the mesons into their constituent quarks. This is an analogue of an insulator- 
metal transition, since the system goes from being an insulator with zero current (in the 
applied field) to a conductor with free charge carriers (the quarks). At finite temperature 
this phenomenon persists, and the dissociation transition becomes subsumed into the more 
familiar meson melting transition. Here, the dissociation phenomenon reduces the critical 
melting temperature. 
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1. Introduction 

In the series of studies of the dynamics of (non-backreacting) quarks [1] in large N SU{N) 
gauge theory using holographically dual[2] geometries in string theory[3, 4, 5, 6], the non- 
perturbative phenomena that can occur in the presence of external electromagnetic fields 
have recently become recognized as readily amenable to computation. The case of external 
magnetic field at zero temperature was studied in refs.[7, 8], where a number of interesting 
phenomena (such as spontaneous symmetry breaking) were readily extracted. There was also 
a recent study of quarks at non-zero baryon number density (see refs.[9, 10]) in the presence 
of a background electric field, where properties (such as the conductivity) of the system were 
uncovered, again using holographic techniques. Very recently, the finite temperature case 
of external magnetic field was studied in some detail, with the uncovering of an interesting 
phase diagram arising from the competing effects of temperature and magnetic field [11]. 

In this paper, we study the same theory (SU{N) gauge theory with non-backreacting 
hypermultiplet quark flavours [1]) — an M = 2 supersymmetric theory in its ground state — 
in the presence of a background electric field, at both zero and finite temperature^. Right at 
the outset, we might expect very different phase structure than the magnetic case since the 
electric field and the finite temperature are working together. The quarks and anti-quarks 
that are bound together to make mesons are oppositely charged under the electric field, and 
so its presence should serve to reduce that binding energy, making it easier for finite tem- 
perature to shake the mesons apart. This is indeed what we see. Moreover, we can quite 

'^We note that another group will present results in this area in a paper to appear shortly [12]. 
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readily identify a sharp first order phase transition even at zero temperature corresponding 
to the dissociation of the mesons into their constituent quarks above a critical value of the 
electric field. Furthermore, since the mesons are electrically neutral, the dissociation repre- 
sents a rather clean example of a "metal-insulator" phase transition, where charge carriers 
are liberated, creating a current in response to the applied field. 

The paper is organised as follows. In section two we introduce the dual string background 
we are going to use throughout the paper and the D7-brane probes that represent the quarks. 
By switching on gauge fields of a particular sort on the world-volume, we can introduce a 
background electric field. This same gauge field may also source a non-trivial current in the 
gauge theory, and we allow for this possibility in our computations. In the next section we 
solve the equations of motion that result from our apparatus, and observe the aforementioned 
phase structure. In particular, when the dissociated or melted phase is present, the system 
has a non-zero result for the response current, representing the flow of the liberated quarks, 
the charge carriers. We map out the phase structure of the system in this section. There 
is an important subtlety that we observe. At high enough values of the electric field, some 
of our solutions for the embeddings of the D7-branes have conical singularities that we can 
characterize, but do not fully understand. It is a possibility that these solutions are corrected 
in the full string theory. However, we believe that their ultimate fate does not change the 
phase diagram we present. Until we know conclusively the fate of these solutions, we are 
careful to clearly show how they arise and where they lie in the phase diagram which we 
uncover. It is to be expected that future work in this area may yield techniques that can 
help further refine our embeddings. In section four we observe some aspects of the meson 
spectrum of the system, and we briefly conclude in section five. 

2. The String Background 

We consider the following form for the metric of the AdSs-SchwarzschildxS*^ background: 

ds'^/a' = — ^ (~Ldt^ + fdxA + ^dr^ + B? cos^ dd^l 
I Ar^E? \ f J 

+R^de'^ + sin^ edfp"^ , (2.1) 

/ = 4r^ - 6\ / = 4r^ + 6^ . 

The AdS5 time and space coordinates are given hy t = and x = respectively 
(which the dual gauge theory will also have) and r G [0,00), which is the AdS radial coordi- 
nate, which we will sometimes denote as x^. We use standard polar coordinates on the S^, 
given by: 

dnl = de^ + cos^ edn'^ + sln^ ^-^-^^ 



2 j„/,2 , „„2<ip , „■ 2/17 ,2 



and dni = d^p^ + cos^ i/jdp + sin^ 61^7^ . (2.2) 



The scale R determines the gauge theory 't Hooft coupling according to R^ = a'^gy-^^^N. 
The parameter b sets the radius of the horizon of the black hole, which in turn sets the 
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temperature, according to T = b/irR^. The coordinate r is related to the traditional radial 
coordinate u by the coordinate transformation r'^ = ^ (^u^ + — 6^^ . To proceed, we 
consider the action for the probe D7-brane to second order in a': 

-irr~ = '-'DBI + '-'WZ 
A// 

= -rD7 I detV2 ^ p [p^^] ^ 27ra'F,b) 

+ (27ray ^ J F(2) A F(2) A P [^4)] , (2-3) 

where P [Gab], P [Bab], and P [C(4)] are the pullback of the background metric, the i?-field, 
and the 4-form potential (sourced by the N D3-branes) onto the D7-brane worldvolume 
respectively, F is the field strength on the D7-brane worldvolume, and Tdy = fJ-j/gs- In 
order to introduce our background electric field (and current), we consider an ansatz for the 
world-volume gauge field of the form [9] : 

Ai{r) = -Et + B{r) . (2.4) 

This ensures a constant electric field Fqi = E \n the gauge theory; B{r) sources a current 
on the world-volume [9]. For the D7-brane embedding, we consider an ansatz of the form: 

e = 9{r) . (2.5) 

The resulting effective action for the D7-brane is: 



Sb7 = -Tj^rNf I "-^^ Pf (1 + r'e'irf) (2.6) 



-4 (27ra')'r4{-/2/S'(r)2 +4/2/^4^2 (1 + ^2^/(^)^ 



1/2 



It is convenient to define dimensionless quantities f, B and E via: 

r = br, B{r) = ^B{f) , E = i^^^E , 0(r) = 0(f) , (2.7) 
which reduces the action to: 

5d7 = -2^ViV/TD76^ [ dt [ dr C , (2.8) 



where C 



cos 9{r) 



16f^ 

with g = Ar'^ - 1 , g = 4f^ + 1 



5^ (g^ - IQf^EA (1 + f'^e'{rf) + Af^g'^gB'{ff 



1/2 



The equations of motion derived from equation (2.8) for B{r) introduces a constant of mo- 
tion T. We can invert the equation of motion to write a first order differential equation for 
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B[r) in terms of T: 



^rgJ (g"^ - I&t-'^eA (1 + f'^e'(fY) 
B'{r) = -f ^ = . (2.9) 

gjg (64f6 (^S'^g'^ - f^g^^ + g'^p cos6 0{f 

In the limit of f ^ oo, the asymptotic solution for B{r) is given by: 

f 

lim B{f) = 6 + — + . . . 
f-»oo 2r^ 

Normalisability of the solution requires us to take 6 = 0. The constant T is related to the 
(vacuum expectation value) vev of the current by the following relation [9]: 

(ji) = {^-f^'^) = -ATr'^^a'b'^VNfTBif , (2.10) 

where V is the spatial volume of the gauge theory. It is convenient to exchange the fields 
B{r) for the constant T by performing a Legendre transformation, defining a new action lur- 

Im = Sj,r - [ fJ-^ 

J Ol<ii 

dC \ 



-2Tr'^V NfT^jb^ J dt J dr - B'{f) 
-27rViV/TD76'^ y dt j df 



v/(l + f2g/(f)2) 
IGr^g^/g 



X y (^52 _ 16f4^2 j |^_64f 6t2^2 + g2g3 ^^^6 

= -27rViV/TD76'lD7 , (2.11) 

where F41 = drAi. It is a simple check to show that SI-dj/ST = B'{f). The expression under 
the square root in equation (2.11) must be positive in order to keep the action real. Either 
both terms should be negative or both should be positive, and they must change sign at the 
same radial distance f^,. This results in two conditions: 

(4f^-l)'-l6r1^2^0 , (2.12) 

- 64f^f 2 (4f^ + 1)^ + [Aft - if [Aft + 1)^ cos*^ 0(n) = . (2.13) 
Equation (2.12) determines in terms of E: 

r, E+ x/PTT 
rl = , (2.14) 
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and equation (2.13) yields T in terms of f^,, and 9{r^,): 

f2 ^ (4r1-l)^(4f4 + l)=^cos6g(f,) 
64f6(4f4 + l)2 

Note that equation (2.15) relates the current to the electric field in a form of Ohm's law, 
and this was used in ref. [9] to determine the conductivity of the fundamental matter in the 
quark-gluon plasma. We refer to f* as the "vanishing locus". From the action l£n defined 
in equation (2.11), we can derive the equation of motion for the embeddings 6{r). Its exact 
form is not illuminating, and so we do not display it here. In the limit of large r, the equation 
of motion asymptotes to: 

{r^9'{r)) + 3r^9{r) = , (2.16) 



dr 



which has as a solution: 



m = ? + (2.17) 



The constants m and c are related to the bare quark mass and the condensate respectively, in 
a manner that is by now very standard [1, 13, 14, 15]. In our notation, the exact relationship 
is given by: 

mg = , (^V') = -Sn'^a'VNjTDrb'^c . (2.18) 



3. Properties of the Solutions 
3.1 Exact Results at Large Mass 

It is instructive to study the properties of the quark condensate as a function of the bare quark 
mass for large mass. This corresponds to m ^ in terms of the dimensionless quantities 
introduced in the previous section. This limit constrains us to consider only the so-called 
"Minkowski" embeddings for the probe D7-brane for which the vev of the current (defined 
in equation (2.10)) vanishes by virtue of equation (2.15). To extract the behaviour of quark 
mass and condensate we linearize the equation of motion obtained from equation (2.11) in the 
same way as described in ref. [11]. Rewriting the result in terms of dimensionful parameters, 
we find the following analytic behaviour for the condensate: 

^^^)--^=^^ ' (3-1) 

where E = lira'E. This may suggest that for high enough bare quark mass, the condensate 
vanishes. However unlike the results for an external magnetic field [11], the condensate 
approaches zero on the positive side of the condensate axis. We verify these results in the 
next section using numerical techniques. 
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3.2 The Case of Vanishing Temperature 

It is instructive to first study the zero temperature case, since a number of important features 
appear in this case. For this simple case, the energy scale is set by -RV^, and therefore it is 
convenient to introduce the dimensionless quantities f and rh via: 

r = RVEr; m = Ry/Em. (3.2) 

We solve the equation of motion for 9{f) using a shooting technique by imposing infrared 
boundary conditions described in refs. [15, 16, 17]. For "Minkowski" embeddings, in order to 
avoid a conical singularity at 9 = 7r/2 [17], we impose the following boundary condition: 

0'{r)\g=./2 = -oo . (3.3) 

For embeddings reaching the pseudo-horizon we find it numerically convenient to shoot for- 
ward (towards infinity) and backward (towards the origin) from the vanishing locus = 1. 
The boundary condition ensuring the smoothness of solutions across the vanishing locus can 
be determined from the equation of motion itself, and it is found to be: 

0'{h) = ^^; eo^ein); (3.4) 
or ou 



de 1 00 
ou u 2 



«(r,) 



From the above boundary condition we can see that the embeddings reaching the vanishing 
locus at 9 = 7r/2 avoids the conical singularity by having a diverging derivative. Next we 
proceed to discuss our numerical findings. 

Numerical analysis shows that there are three different classes of embeddings, which can 
be classified by the topology of the probe D7-brane. First, we have the smooth "Minkowski" 
embeddings that close above the vanishing locus, due to the shrinking of the S"^ that the 
probe brane wraps. Second, we have the embeddings that reach the vanishing locus before 
the shrinks. These embeddings can in turn be classified into two different categories. 
The first are singular solutions that cross the vanishing locus and close before reaching the 
origin. These have a conical singularity at the closing point because the derivative does not 
diverge. We discuss these solutions further below. The second are smooth solutions that pass 
through the vanishing locus and reach the origin with no singular behaviour. These different 
embeddings are summarized in figures 1 and 2. We can see explicitly the conical singularity 
appearing for the brane (second curve in figure 2) since it intersects the vertical axis before 
reaching the event horizon. This class of singular solutions starts appearing after a certain 
value of 6*0 = (^'o)min t^^t depends on the magnitude of the electric field; this value is close 
to 7r/2 (which corresponds to the critical embedding), and the singular solution persists until 
00 = 7r/2. 

The existence of these embeddings suggests the possibility of a transition in topology 
of the probe brane as a function of the bare quark mass, as is well known to exist at finite 
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Figure 1: The dashed red curves represent solutions smoothly closing before reaching the vanishing 
locus. The solid blue remaining curves correspond to embeddings passing the vanishing locus (denoted 
by a semi-circular dashed black curve). 




Figure 2: Two kinds of embeddings that pass the vanishing locus corresponding to D7-brane closing 
at the origin (left), and forming a conical singularity (right) respectively. The semi -circular dashed 
(black) curve corresponds to the vanishing locus. 

temperature. Here, the external electric field alone induces the transition. The Minkowski 
embeddings simply have a shrinking 5^, while the embeddings reaching the origin are distin- 
guished by having an S"^ shrinking as well as touching the AdS horizon. However the presence 
of singular solutions makes the nature of the transition subtle, as we discuss later. We later 
see that the finite temperature case also shares a similar classification of embeddings. 

From the asymptotic behaviour of these embeddings we can extract the condensate as 
a function of the bare quark mass. In figure 3 we show this dependence. There are two 
important mass scales. One is where the phase transition between the two types of embedding 
occurs, which we label rhcr, and the other is where the singular solutions first appear, which 
we label m* . These values are in turn determined by the parameters (0o)cr and (0o)min, 
respectively. 

In figure 3, the multi-valuedness is indicative of a phase transition, as was speculated 
earlier based on the general arguments. Furthermore we observe that there is no chiral 
symmetry breaking, since in all cases the condensate vanishes for vanishing mass. This is 
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Figure 3: The condensate as a function of bare quark mass to. The curve segment coming in 
from the left (sohd blue) corresponds to embeddings reaching the vanishing locus, and the curve 
segments going out to the right (dashed red) correspond to Minkowski embeddings. There is a family 
of solutions that contain a conical singularity between the vanishing locus and the origin. They exist 
for (0o)min < ^0 < (^o)cr- On the right is a magnification of the turn-around region where these 
segments join, showing multi-valuedness. Since it has been a good guide in earlier work [15], we use 
an equal-area law argument to determine the presence of a first order phase transition at TOcr where 
there is a jump from one type of embedding to another. 



to be contrasted with the case of having an external magnetic field where chiral symmetry 
breaking is observed [11]. 

The latter observation fits our intuition that the mesons in this theory have a binding 
energy that is proportional to the constituent quark mass [13]. For a given quark mass, 
the introduction of an electric field reduces the binding energy of the quarks. This causes 
two things to happen: First, it inhibits the formation of a chiral condensate, and second, it 
ultimately dissociates the mesons into its constituent quarks. This dissociation is in fact a 
transition from an insulating to a conducting phase, mediated by the external electric field. 
On the dual gravity side this corresponds to a transition from Minkowski embeddings to 
embeddings that reach the origin. 

We must note that the appearance of the singular solutions (those that have a conical 
singularity) are not well understood at the moment. There is therefore the possibility that 
there is as yet not identified intermediate phase right after the mesons dissociate. We can 
make our dissociation transition explicit. From equation (3.2) we can deduce the exact 
dependence of m* and mcr on the electric field E as m* = R\f^m* and mcr = R\f^ rhcr- 
The latter is the value of the critical mass for a given electric field. Using an equal-area law 
(see ref. [15])^, we can also independently (as a test of our methods for later) determine this 
dependence numerically for our solutions, and we display this in figure 4. We also include 

^It would be interesting to explore comparing the results of the equal-area prescription with the results of 
the thermodynamic potential computation using an appropriate brane probe action (DBI supplemented with 
terms corresponding to the response to the external fields.). This we will leave to a future project. 
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the critical mass for which the singular solutions appear as a function of the field. We see 
that the analytic behaviour deduced above is nicely confirmed by the numerics. Indeed the 
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Figure 4: The phase diagram showing the electric field induced first order phase transition separating 
the dissociated phase from the phase with stable mesons. The blue (solid) curve separates the phases. 
The region within the blue (solid) curve and red (dashed) curve represents the region where the 
embeddings with a conical singularity become thermodynamically relevant. Note that rrici- has units 
of R and E is dimensionless. 

singular solutions seem to be unavoidable in the phase diagram (at least at the level of our 
analysis). However they seem to lie very close to the dissociation transition and, in the limits 
S — > and E oo, the two curves tend to merge. 

Presently we do not completely understand the role played by the singular solutions in 
the physics. One simple possibility (but not the only one) is that stringy corrections smooth 
out the singularity in the interior, while preserving the asymptotic behaviour. This would 
mean that we have a simpler phase diagram than that in figure 4 (i.e., with no dashed line). 
In fact, a study of the meson spectrum appears to hint at the fact that the gauge theory is 
actually "blind" to the physics below the vanishing locus, at least for the type of questions 
we are asking. We elaborate on this in section 4. 

We next turn to the finite temperature case. It is worth noting that a simple check of 
our results in the next section is that the large electric field (compared to temperature) limit, 
should recover the physics that we have seen here, and indeed it does. 

3.3 The Case of Finite Temperature 

In order to proceed, we again solve the equation of motion for 9(r) numerically using Mathe- 
matica. We use a similar shooting technique and boundary conditions outlined in the previous 
section. For Minkowski embeddings we impose equation (3.3). For black hole embeddings, 
the appropriate initial condition to ensure smoothness of the embedding when we stitch the 
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solutions across r* is: 



du do 
df du 



(3.5) 



86 _ (3m4 - 1) - V(3m^ - 1)2 + 9u^{u^ - 1) tan2 



3{i(M4 - 1) tan 6*0 

We show several solutions for £^ = 1 in figure 5. We can extract from these embeddings the 
r sin 6 



1.5. 
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Figure 5: Several solutions for classical D7-brane embeddings with E = 1. The inner semi-circle 
corresponds to the event horizon, and the outer semi-circle corresponds to the vanishing locus. 

condensate c as a function of the bare quark mass m by analyzing the asymptotic behavior 
of the solutions. We show several of these solutions for various dimensionless electric field 
values in figure 6. The first order phase transition described in refs. [14, 15] persists in the 
presence of an electric field. We determine the critical mass at which this transition occurs 
using an equal-area law described in ref. [15], and we present the phase diagram in figure 7. 
The numerical results confirm our intuition: at zero temperature, the mass of the meson is 
proportional to the bare quark mass, and hence the binding energy is proportional to the bare 
quark mass [13]. Turning on an electric field decreases the binding energy (the electric field 
will pull the quark-antiquark pairs apart), and only mesons built out of sufficiently heavy 
quarks-i.e. above a critical mass-will survive. At finite temperature, we can think of the 
increasing of the critical mass as equivalent to the decreasing of the critical temperature. 
Therefore, lowering the binding energy causes mesons to melt at lower critical temperatures. 
In addition, since c = c/6^ the magnitude of the dimensionless condensate increases with 
decreasing critical temperature (as we see in figure 6). Furthermore, the stronger the electric 
field, the more initial binding energy is required and therefore a higher critical mass is needed. 

We can deduce the asymptotic behavior of the phase diagram by considering large electric 
fields. The dominant energy scale is set by the electric field, and therefore in that limit we 
are effectively at zero temperature. To see this, one should note that the energy scale in the 
geometry is set by r^,, and for strong electric field, r^, is much larger than the position of the 
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Figure 6: The condensate as a function of bare quark mass for increasing electric field. From bottom 
to top, the electric field values are E = 0.01, 0.1, 0.3, 0.5, 0.7, 0.9. 

event horizon. Therefore, we should expect to reproduce the results for pure AdSs x in 
an external electric field described in the previous section. In the case of pure AdSs x in 
an external electric field, the only energy scale is given by R^/E, and therefore we can easily 
show that the bare quark mass and condensate scale as: 

m oc RsfE , 

c oc R^E^/^ . (3.6) 

This result suggests that: 

l/mcrit « • (3.7) 

This dependence was checked numerically (and shown in figure 7) and appear to be valid 
for E > 2. It is instructive to compare the phase transition to the one for zero electric 
field and finite temperature. In the latter case the Minkowski and black hole embeddings are 
separated by a "critical" embedding with a conical singularity at the event horizon as shown 
in figure 8 . The straightforward generalization of this picture to the case of finite electric field 
suggests the existence of a "critical" embedding with a conical singularity at the vanishing 
locus (r = n). However, similarly to the zero temperature case, a more detailed study of 
the transition from Minkowski to black hole embeddings reveals an interesting third class of 
embeddings. These are embeddings which enter the vanishing locus at r^,, and have a conical 
singularity above the event horizon. See figure 9. These solutions correspond to embeddings 
which enter the vanishing locus at close to 7r/2. Such solutions extend from 6q = 7r/2 
to some minimal value (^o)min) below which the embeddings have no conical singularity and 
simply fall into the black hole. In figure 11 one can see examples of embedding for different 
values of ^o- Finite temperature adds richness to the picture by adding the possibility that 
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1 2 3 4 5 

Figure 7: The phase diagram depicting the first order phase transition, denoted by a heavy line 
(blue) that separates the melted and stable meson phase. The thin black curve is the best fit curve 
given by const, x E^^/"^ for large electric field. 




Figure 8: For zero electric field and finite temperature, the Minkowski and black hole solutions are 
separated by a "critical" embedding, which has conical singularity at the event horizon. 



these singular solutions may or may not be bypassed by the otherwise melting transition 
induced by temperature. This depends entirely on the value of the electric field for a given 
temperature, and for sufficiently high electric field the singular solutions appears unavoidable 
(at least at this level of analysis — as mentioned earlier, there may well be stringy corrections 
to the conical geometry in the interior). Summarizing the two possible scenarios mentioned 
above, we have: 

1. For sufficiently weak electric field or large temperature {E < Ecr) the {9o)mm embedding 
is in the vicinity of the "critical" (^o = ^/2) embedding and is bypassed by the phase 
transition. The thermodynamically stable phases correspond to smooth Minkowski and 
black hole embeddings on the gravity side or meson gas and quark gluon plasma states 
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Figure 9: For a finite electric field there are two "critical" solutions corresponding to embeddings 
with a conical singularity at the vanishing locus (r = r*) and at the event horizon respectively. The 
intermediate solutions correspond to black hole embeddings, which have a conical singularity between 
the event horizon and the vanishing locus. 
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Figure 10: The physical scenarios described in the text for (left) E < and (right) E > Ea-- To 
the left, we have (^o)cr ^ (^o)inin, so the phase transition bypasses the singular branch of solutions. To 
the right, we have (0o)cr ^ (^o)mini so the singular solutions appear on the resulting phase diagram. 
Further discussion is available in the text. 



in the dual gauge theory. 

2. For E > Ecr the (6'o) min embedding is thermodynamically stable and corresponds to 
some bare quark mass rh* < rhcr- This means that some of the classically stable 
black hole solutions, namely the one corresponding to bare quark mass in the range 
fh* < m < rhcr, will have a conical singularity. 

There is an example of each such scenario in figure 10. This suggests a richer structure 
for the phase diagram than the one presented in figure 7 and in particular the existence of 
a special point at E = E^r and ?ficr(-E'cr)- The corresponding phase diagram is presented in 
figure 12. The new choice of dimensionless variables on the axes is made in order to emphasize 
the existence of the non-smooth {i.e., conically singular) black hole embeddings. The solid 
curve corresponds to the phase curve from figure 7 (but note that the axes are different here, 
and hence the shape) and the dashed curve separates embeddings with a conical singularity 
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Figure 11: Examples of D7-brane embeddings for > (^o)min, do = (^'o)min and 6*0 < (0o)min 
respectively. The inner semi-circular segments represent the event horizon, while the outer semi- 
circular arcs correspond to the vanishing locus at r = . The numerical value of the minimal angle is 
(^o)min ~ 1-345 and the parameter E — 3. 



from smooth black hole embeddings (that lie below it). The area above the solid curve is 
in the stable meson phase and the area below the dashed curve is in the melted/dissociated 
phase. The area between the curves corresponds to the embeddings that have a conical 
singularity before ending on the horizon. The vertical dashed line corresponds to the critical 
value Eci ~ 1.26, below which the conical solutions are not part of the energetically favourable 
solutions and so do not complicate the story. In order to relate the results in figure 12 to 
those of the zero temperature case, we first note that the T = dimensionless variable rh can 
be related to the finite T variable rh via: 

m m 

m = j= = lim — = . (3.8) 

rVe t^o ^/^ 



In addition, if ^ = b'^E/R^ is to remain finite as 6 — > 0, must go to infinity. Therefore, 
the T ^ limit corresponds to studying E ^ co in figure 12. We see that the phase curves 
from figure 12 qualitatively approach the zero temperature (large electric field) case of the 
previous section, with the critical parameters m* « 1.22 and rhcr ~ 1.32. 



4. Meson Spectroscopy 

We end by making some observations about the meson spectrum. To proceed in studying 
the meson spectrum, we first absorb our gauge field into the B field to emphasize that it is 
of the same order in a' at the classical level of our Dirac-Born-Infeld action. This simplifies 
our notation since now denotes the second order in a' gauge field. To study the meson 
spectrum, we need to consider quadratic (in a') fluctuations about our classical embeddings 
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Figure 12: The solid curve correspond to the phase curve from figure 7 and the dashed curve separates 
embcddings with a conical singularity from smooth black-hole embeddings. The area above the solid 
curve is in the meson gas phase and the area below the dashed curve is in the quark-gluon plasma 
phase. The area between the sohd and dashed curves corresponds to embeddings which reach the black 
hole horizon but have a conical singularity before reaching the event horizon. Their physical status 
remains to be determined but may be irrelevant for this particular phase diagram (see text). The 
vertical dashed line corresponds to the critical value E^r ~ 1.26 where the conical solutions become 
thermodynamically relevant. 



for the D7-branes [1]. To this goal, the relevant pieces of the action are: 

5 = -— /" d^iyJgah + Bab + 2T:a'Fab + (2W) ^l^ [ F(2) A B^^) A P 
9s J J Ms, 



where 



+ (27ra')V7^/ F(2) A F(2) A P [C( 
^ J Ms 



(4) 



(4) 



'Ms 
1 

9s 



4)J , 



dt A dx^ A dx"^ A dx^ , 

(l — cos^ 9) sin i/j cos il^ dil) f\ da f\ d[) /\ • 



(4) 



(4.1) 



To study the quadratic fluctuations, we consider an expansion for the transverse fields of the 
form: 



eo(r) + 27ra'x(r) 
+ 2W$(^") , 



(4.2) 
(4.3) 
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where is the classical embedding found in previous sections. Expanding the action in 
equation (4.1) to second order in 27ra', we find the relevant terms for our fluctuation analysis 
are: 

1 I — fj^^ „ „ 1 



— Cp2 



^, -EE^'GeedaXObX ^ 
1 



+ 2^" 

1 r 

2 

1 



dl^TE - du (e'GeeE'^^deV^E ' 
EG^^E'^^da^db^ , 



-ES'^'S^'^FbcFda 



ab 



R'du 

3 



-E 



C 



wz 

F-F 



-4R cos 9 sin 9 sin ip cos ^/;i?oi-^23'^ 

„.4 



7? VP c 



mnop 



(4.4) 
(4.5) 
(4.6) 

(4.7) 
(4.8) 
(4.9) 



Here Eat = Qab + Bab-, where gab is the induced metric at the classical level. G will be reserved 
for the space-time metric. In addition, we use that E'^^ = S""^ + J"'', where 3°"^ = S^"" and 
jah ^ _jba_ rp^g indices a,6,c,d range over the world volume coordinates, and the indices 
m,n,o,p range over the directions 4,5,6,7. We have reverted to using the coordinate u 
instead of r since the equations are simpler in this coordinate. The equations of motion 
derived from these Lagrangian terms allow us to consistently take Ai{t,u) ^ 0, A4^{t,u) ^ 0, 
and the remaining gauge fields to be zero. This implies that the field x is coupled to Ai 
and A4, whereas the field ^ decouples from the gauge fields. Let us begin by studying the 
spectrum for the zero temperature case. For simplicity, we consider fluctuations about the 
9 = embedding, in which case, the field x decouples from the gauge field. We consider an 
ansatz of the form: 

X = x(^i,t)=e(n)e-^^S (4.10) 



where u = uRy E, t = ty E/R, and 6j = ujR/y E are dimensionless quantities. We would 
like to focus our attention to the equation of motion near the vanishing locus, in this case, at 
u = = 1. The equation of motion has the form: 



e"(n) + 



1 



l + i 



V6J 



+ 



1 



u 



Co 5u) 
\/6 



12 



. 



(4.11) 



We emphasize that this equation assumes that T is positive. We find that this equation 
satisfies an ingoing wave solution of the form: 

C(n) = (n-n,)-*^/^C(n) , (4.12) 

The outgoing wave solution corresponds to simply taking the complex conjugate of the ingoing 
wave solution: 

Xout = Xin • (4-13) 
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In examples involving string worldsheets [18, 19], such behavior was interpreted as an effective 
horizon for the worldsheet, and only the ingoing wave solution was considered physical. Here, 
the ingoing solution implies that embeddings that touch the vanishing locus must fall towards 
the event horizon. This is simply the statement that the conducting phase can only be 
associated with black hole embeddings. We emphasize that the situation here is not exactly 
the same as in ref. [18, 19]; the vanishing locus does not correspond to the horizon of the 
induced D7-brane metric. However, there are indications that the gauge theory only sees the 
physics up to the vanishing locus. The current T is defined at n*; the asymptotic behavior 
of the branes can be solved for by imposing initial conditions at u*, so the bare quark mass 
and condensate can be identified without having to solve for the embedding up to the event 
horizon. From our results above, the same holds true for the meson spectrum. Therefore, one 
may be tempted to interpret the vanishing locus as an effective horizon; however, Ramond- 
Ramond (RR) charge conservation on the D7-brane requires us to continue the embedding 
up to the physical event horizon. Perhaps for the questions we are asking, the gauge theory 
is ignorant of the physics inside the vanishing locus, and hence the conical singularities we 
have encountered and their resolution do not alter the phase diagrams we have presented. 

Therefore, let us proceed and check the behavior of the fluctuations for the T 7^ case. 
For a similar ansatz as before, near the vanishing locus, we again find ingoing wave solutions 
of the form: 

e(n) = (n-n,)-^'^''C(n) , (4.14) 

where u = u/h^ i = bt/R"^, Co = R'^oj/b, and a = ^ 2'U^^{2+3^^f ' ^* "i^^ermines the position of 

the vanishing locus and is given by \/l + E'^. Therefore, we find that the solution becomes 
exactly the incoming-wave quasinormal mode solution studied elsewhere [20, 21, 22, 15]. 
Therefore, this result confirms that we still have a "melted" phase for the mesons correspond- 
ing to black hole embeddings. 



5. Conclusions 



It is very encouraging that non-trivial non-perturbative phenomena resulting from exter- 
nal fields such as those we have seen here (a dissociation phase transition, metal-insulator 
transition and the associated response current) can be so readily extracted in this kind of 
holographic study. We found that since the electric field works together with the presence 
of finite temperature, the resulting phase diagram which accounts for the effects of both is 
rather simple (ignoring the complication of the special conical solutions — see below). This 
paper is a natural companion to our recent work [11] on the same system in an external mag- 
netic field, where we mapped out somewhat richer phase structure, seeing how magnetically 
induced spontaneous chiral symmetry breaking (studied in this context in ref. [7]) and the 
meson melting fit together in the diagram. 

Clearly, the story is not quite complete in this electric case, since above a certain value of 
the electric field, some of our solutions develop conical singularities in the interior. We would 
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like to know more about the nature of these solutions. As already stated, one possibility 
is that the solutions are locally (near the singularity) corrected by stringy physics, perhaps 
smoothing the conical points into throats that connect to the horizon. We would expect in 
that case that our phase diagram would largely remain intact, since the values of the masses 
and condensates for each solution are read off at infinity, far from where the conical singularity 
develops. The results for the relative free energies of those solution (compared to the other 
solutions at the same values of the mass) would be the same, and so the complete story would 
be unaffected. Furthermore as noted in section 4, key aspects of the meson spectrum lend 
support to the idea that the phase diagram is insensitive to the local details of the conically 
singular features. It would be of value to find methods for making this indication more robust. 
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